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AGN heating in the centres of galaxy groups: a statistical study 
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ABSTRACT 

We present gas temperature, density, entropy and cooling time profiles for the cores of a 
sample of 15 galaxy groups observed with Chandra. We find that the entropy profiles follow a 
power-law profile down to very small fractions of i?5oo- Differences between the gas profiles 
of groups with radio loud and radio quiet BGGs are only marginally significant, and there 
is only a small difference in the Lx ■ Tx relations, for the central regions we study with 
Chandra, between the radio-loud and radio-quiet objects in our sample, in contrast to the 
much larger difference found on scales of the whole group in earlier work. However, there is 
evidence, from splitting the sample based on the mass of the central black holes, that repeated 
outbursts of AGN activity may have a long term cumulative effect on the entropy profiles. 
We argue that, to first-order, energy injection from radio sources does not change the global 
structure of the gas in the cores of groups, although it can displace gas on a local level. In 
most systems, it appears that AGN energy injection serves primarily to counter the effects 
of radiative cooling, rather than being responsible for the similarity breaking between groups 
and clusters. 
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the inflated bubble will rise due to buoyancy, displacing gas fur- 
ther. It is thought that if repeated cycles of radio outbursts occur, 
and the effects are integrated over time, then radio source heating 
may co unteract catastro phic cooling, even when no active source is 
present teegelmaj200lh . 

It has also been suggested that the energy input into the 
IGM by a radio-loud AGN may also be responsible for similar- 



1 INTRODUCTION 



Over the past decade cosmological simulations have become in- 
creasingly sophisticated to the point where we can now follow the 
forma tion and evolution of individual galaxy groups and clusters 
(e.g. iThomas et al1l200lh . Our knowledge about the evolution of 
groups and clusters has also increased, due to a combination of de- 
tailed X-ray studies using ROSAT, Chandra and XMM-Newton and 
deep, low fre quency radio studie s using the VLA radio telescope 
(for e xample iNulsen et al.l 120051 iHeinz et al.l |2002| . iFabian et al.l 
l200d) . 

However, our understanding of galaxy groups and clusters is 
currently impeded by two 'heating problems' . Whilst recent work 
has shown that the dark matter distribution o f galaxy cl usters scales 
self-s imilarly (e.g. iPointecouteau et al.l 120051 and IVikhlinin et al.l 
2006), the baryonic component (the hot gas) does not follow s elf- 
similar scaling mode ls (see for example. ISanderson et al ]|2003l and 
IPopesso et al. 120051) . This appears to re sult from excess entrop y 
in low mass clusters jPonman et al.ll 199^ and lPonman et al.ll2003l) . 
which might be a consequence of feedback from supernova explo- 
sions or active galactic nuclei (AGN). Secondly, the hot gas in clus- 
ter cores is not coolin g at the rates expecte d from the short cool- 
ing times inferred (e.g. lPeterson et al1 l200ll). suggesting that some 
heating process must act to offset cooling. Since massive early-type 
galaxies are found in the cores of almost all X-ray bright groups and 
clusters, and are believed to host massive black holes, AGN heat- 



ity breaking in group s a nd clusters (for example iKaiserl \\ 99 1[ 
IValageas & SiUj 19991 and lNath & Rovchowdhury||2002l) . Occam's 
razor suggests that both similarity breaking and the prevention of 
catastrophic cooling may be related, since both occur in the same 
systems, and a similar source of energy injection is thought to cause 
both. The best place to look for evidence of such a connection is 
probably not in rich clusters, but galaxy groups, since the gas in 
the shallower potential wells of groups should be more strongly af- 
fected by any embedded heat source, and it is in groups that excess 
entropy is most clearly apparent. 

The recent study of ICroston et al.l d2005l) found that the bulk 
X-ray properties of galaxy groups are related to the presence of ra- 
dio galaxies within them; groups with radio-loud AGN tend to fall 
below the Lx '■ Tx relation, establishing a prima facie case for 
a connection between AGN heating and group scaling properties. 
However, since AGN radio outbursts are thought to be periodic, 
and the way in which the intergalactic gas responds to the out- 
bursts is poorly understood, it is impossible to say whether there are 
long-lasting and cumulative effects on the baryonic components of 
groups (such as would be required to account for similarity break- 
ing) or whether AGN activity simply induces short-lived excursions 
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in the Lx '■ Tx plane with no lasting effects on the temperature 
and entropy of the systems. Further, the inference that it is the ra- 
dio source that has an effect on the IGM, rather than the conditions 
in the IGM being responsible for triggering a radio source, may be 
incorrect; it could be that the state of the IGM and the conditions re- 
quired to trigger a radio source form a feedback loop, such that the 
IGM attains a specific state, which triggers a radio outburst, which 
in turn affects the IGM, and eventually quenches the outburst, al- 
lowing the IGM to return back to the initial state which triggered 
the outburst. 

There appear to be three possible interpretations: 

(i) The radio source could modify the observed X-ray proper- 
ties, through heating or displacing gas. 

(ii) Alternatively, the observed X-ray properties could be caus- 
ing the radio source by providing conditions that favour the trig- 
gering of a radio outburst. For example, the hot g as could provide a 
suitable reservoir of gas for Bondi accretion (e.g. lAllen et al .l2006h . 
or the hot gas may have a steep pressure gradient required to col- 
limate larg e scale jets such as t hose found in 3C31, in the group 
NGC 383 <Laing & Bridlell2002h . 

(iii) There is a third property of the system, such as BGG mass 
(see Section [3~5l >, or overall mass of the group, which correlates 
with both the radio and X-ray properties, and is responsible for 
both. 

In this paper, we use a sample of 15 galaxy groups to search 
for evidence of radio source and hot gas interaction at the centres 
of galaxy groups and to try to distinguish between the interpre- 
tations di scussed above. The groups were taken from the GEMS 
sample of lOsmond & Ponmanl ( I2004I) . and those groups in their G 
and H samples that had good quality archival Chandra data avail- 
able were used for this study (and correspond to a subset of the 
ICroston et al.l d2005h sample). It should be noted that the parent 
GEMS sample is not, in itself, unbiased, favouring systems with 
higher X-ray luminosities. In this study, we have used a sub- 
set of the GEMS sample which had available archival Chandra 
data with observations longer than 30 ks; thus we do not expect 
to explicitly bias our sample further. However, from Fig. [TJ it 
can be seen that compared to the sample of lCroston et alj 1 2005) 
our selection procedure does tend to favour higher luminosity 
systems. Since the X-ray luminosity range in which the ROSAT 
sample shows effects of heating is nevertheless adequately sam- 
pled by our Chandra sample, this slight bias should not affect 
our conclusions. 

We present our sample and outline our data analysis in Sec- 
tion |2j we discuss our results in Section [3] and present our con- 
clusions in Section [5] We assume that Ho — 72 km s _1 Mpc - 
throughout. 



2 DATA ANALYSIS 

2.1 Details of the sample and data reduction 

Details of the groups used in this study are given in Table Q] Each 
Chandra dataset was cleaned following the procedure outlined in 
the CIAO threads on-line; the data were reprocessed to apply the lat- 
est calibration products, and a new events file was created. Periods 
of high background were then removed from the data by creating a 
light-curve of the background, and using a 3-a clipping algorithm 
to identify and remove times of high background. These datasets 
were then used in the subsequent spectral and imaging analysis. 
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Figure 1. The ROSAT Lx '■ Tx relation of the groups used in the Croston 
et al (2005) study. The points marked with stars indicate the groups analysed 
in this study; whilst higher luminosity systems are favoured, this should not 
effect our conclusions (see text). 



2.2 Imaging 

For each group, the cleaned data were smoothed using the adap- 
tive smoothing algorithm CSMOOTH which scales the size of the 
smoothing Gaussian function until a pre-defined number of counts, 
which corresponds to a user specified significance, is met. This 
implies that only structure which is 'real' in a statistical sense is 
smoothed. 

From the images, there does appear to be evidence of dis- 
turbed hot gas in so me group cores; g roups such as NGC 741 
(Fig.|2j and HCG 62 jBirzan et al]|2004?) show 'ghost' X-ray cav- 
ities, and NGC 4261 shows 'filaments' of hot gas outlining the ra- 
dio lobes (shown in Fig.|2j, which also extend to larger scales (see 
ICroston et al]|2005h . If the hot gas shows evidence of spatial inter- 
action with the radio plasma injected into the group by the AGN, 
then it may be the case that a change in radial temperature, density, 
or entropy distributions could be seen if the AGN is injecting extra 
energy into the group. 



2.3 Spectral analysis 

2.3.1 Obtaining temperatures and luminosities 

In order to investigate the effect of radio outbursts on the Lx '-Tx 
relation for the inner regions of galaxy groups, spectra were ex- 
tracted for the systems in the sample out to a fiducial radius of 
0.05i?5oo- This radius was used as we had data for all groups 
in our sample, apart from NGC 4261 which was excluded from 
the Lx'.Tx analysis since the original X-ray observation did 
not have data out to O.O5.R500. Our measu rement of fl s pp wa s 
obtained using the Rtmq ■ Tx relation of Willis et al] J2005I) . 
who use the XMM-LSS for a sample of galaxy groups and clus- 
ters. In addition to the groups shown in Table Q] we also obtained 
temperatures and luminosities for a further 4 groups (NGC 1587, 
NGC 3665, NGC 3923, and IC 1459) that were not included in the 
original sample due to the quality of the data (and hence did not 
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Table 1. Details of the groups in this study, and the corresponding Chandra data. 



source name z group centre co-ordinates R500 Chandra obs date observed cleaned exposure 
02000 #2000 (kpc) number time (ks) 



NGC 383 


0.017 


01 


07 


25 


32 24 46 


510 


2147 


2000 


1 1 


06 


42 


8 


NGC 533 


0.019 


01 


25 


31 


01 45 34 


410 


2880 


2002 


07 


28 


37 


6 


NGC 720 


0.0058 


01 


53 


01 


-13 44 20 


260 


492 


2000 


10 


12 


39 


6 


NGC 741 


0.019 


01 


56 


21 


05 37 45 


440 


2223 


2001 


01 


28 


30 


3 


NGC 1407 


0.0059 


03 


40 


12 


-18 34 49 


400 


791 


2000 


08 


16 


48 


6 


NGC 3607 


0.0032 


11 


16 


55 


18 03 04 


200 


2073 


2001 


06 


12 


38 


5 


NGC 4073 


0.020 


12 


04 


27 


01 53 35 


510 


3234 


2002 


11 


24 


30 





NGC 4261 


0.0075 


12 


19 


23 


05 49 29 


460 


834 


2000 


05 


06 


34 


4 


NGC 4325 


0.026 


12 


23 


07 


10 37 15 


350 


3232 


2003 


02 


04 


30 


1 


NGC 4636 


0.0031 


12 


42 


50 


02 41 15 


350 


4415 


2003 


02 


15 


74 


3 


HCG 62 


0.014 


12 


53 


06 


-09 12 14 


490 


921 


2000 


01 


25 


47 


4 


NGC 5044 


0.0090 


13 


15 


24 


-16 23 09 


440 


3225 


2002 


06 


07 


83 


1 


NGC 5171 


0.023 


13 


29 


22 


11 44 06 


410 


3216 


2002 


12 


10 


34 


7 


NGC 5846 


0.0057 


15 


06 


29 


01 36 21 


320 


788 


2000 


05 


24 


29 


9 


NGC 6338 


0.027 


17 


15 


23 


57 24 41 


620 


4194 


2003 


09 


17 


47 


3 



NOTES: R500 is c alculated using the ROSAT temperatures of Osmond and Ponman (2004), and the scaling relation of R500 = 391 X T 0, 
Iwillis et al.l2005l) . 



Table 2. Radio powers, galactic velocity dispersions, <r, and K-band magnitudes, Mj{ for the sources in our sample 



Source 


1 .4 GHz radio power 


radio flux 


size of radio source 


(7 


(7 


M K 




WHz" 1 


reference 


at 1.4GHz(kpc) 


km s — 1 


reference 


mag. 


NGC 383 


3.2 x 10 24 


1 


360 


288.0 


3 


-23.57 


NGC 533 


2.2 x 10 22 


1 


20 


224.0 


4 


-25.84 


NGC 720 


< 1.0 x 10 20 


2 




273.0 


3 


-24.57 


NGC 741 


7.9 x 10 23 


1 


150 


270.0 


5 


-26.13 


NGC 1407 


6.6 x 10 21 


1 


to 


259.7 


6 


-25.15 


NGC 3607 


1.2 x 10 20 


1 


4 


240.0 


7 


-23.49 


NGC 4073 


< 3.9 x 10 20 


2 




276.0 


8 


-26.04 


NGC 4261 


2.3 x 10 24 


1 


80 


316.0 


9 


-25.10 


NGC 4325 


< 6.3 x 10 20 


2 








-24.77 


NGC 4636 


6.1 x 10 21 


1 


3 


208.0 


9 


-23.97 


HCG 62 


3.2 x 10 20 


2 


15 


251.0 


10 


-25.09 


NGC 5044 


6.2 x 10 21 


1 


15 


240.0 


11 


-25.08 


NGC 5171 


< 2.6 x 10 21 


2 








-24.73 


NGC 5846 


1.5 x 10 21 


2 


1 


261.0 


9 


-24.93 


NGC 6338 


9.5 x 10 22 


2 


6 


347.0 


12 


-25.84 



NOTES: Column 3 gives the reference for the radio flux density: 

1 Radio flux from NED. 

2 Radio flux from NVSS fcondon et al.ll998l) . 

Column 6 gives the reference for galactic velocity dispersion: 
SimthetalJfeOOOl) 
Mtiller et aljjl999h 
Bonfantiet al l Q995I) 



Beuing et alj 120021) 



Proctor & Sansodfc002l) 



Fisher et al. 1 199: 



m,: 

>5) 



Bernardi et alj J2002I) 



9 
10 
11 
12 

K-band magnitudes collated from 2MASS 



Ramnazzo et al. 11998 
Carollo etal|h993 
Wegne ret alii 1999 



have a deprojection analysis done, see Section [2.3.2t . Spectra were 
extracted using the CIAO acisspec script, excluding point sources 
but not the central AGN where one was present, and using a local 
background immediately outside the region from where the spec- 
trum was extracted. The spectra were fitted in XSPEC with an ab- 



sorbed MeKaL model, and a power law for the AGN component 
(if present) which was left free to fit the data, and a power law with 
a fixed index of 1.72 to account for any contribution from X-ray bi- 
naries. A temperature and unabsorbed luminosity within 0.05i?5oo 
was measured for each group. The results are shown in Fig[3] 
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Figure 2. Groups in our sample which exhibit spatial evidence of interaction between the hot IGM and the AGN. NGC 4261 shows filaments of hot gas 
outlining the radio lobes, whereas NGC 741 and HCG 62 exhibit ghost cavities, where there is no clear correspondence between AGN emission (if any) and 
features in the X-ray. 



2.3.2 Deprojection 

If the radio source in a specific group is having an effect on the 
bulk properties of the hot gas, rather than simply displacing gas 
and causing local perturbations, then it might be expected that an 
effect would be seen in the overall radial properties of the group. 
If radio source heating is occurring, then elevated temperatures and 
excess entropy may be present in those groups with current activity, 
compared to those which are currently radio quiet. By examining 
the radial profiles of the gas properties, and comparing the profiles 
of radio-loud groups with those of radio quiet groups, differences 
in temperature and entropy may be observed. 

For these reasons, we analysed the galaxy groups using an 
onion skin deprojection method to obtain radially deprojected tem- 
perature profiles. In doing the onion-skin deprojection, we assumed 
that the group was spherically symmetric, and that the gas had a 
filling factor of unity. Spectra were extracted in concentric annuli 
centred on the brightest group galaxy (BGG), with the outermost 
annulus being used as a local background. The width of each annu- 
lus was determined by the number of counts in each annulus; to do 
our spectral analysis, a minimum of 1000 counts was required in 
each annulus. We then followed the deprojection method described 



in IJetha et all J2005h . fitting each spectrum with a MeKaL model 
to obtain temperatures and normalizations for each shell. 

The normalization of the MeKaL model (N me k) is related to 
the density of the gas in the spherical shell by 

-14 



N„ 



10" 



4tt [D a (1 + z)Y 



(1) 



where n e and n p are the electron and proton densities of the gas and 
Da is the angular size distance to the source. Further assuming that 
n p = 1.18n e , and that the gas to be fitted occupies a volume V, 
then the density of the gas is 



4tt [D a (1 + z)] 2 N mek 
1.18 x 10~ 14 V 



(2) 



The density and temperature profiles obtained from the depro- 
jection were then used to calculate entropy, pressure and cooling 
time profiles for the sources. We calculate the entropy index, S(r) 
for a radial bin using: 



S(r) = T(r)n(ry 2/ \ 



(3) 



where T is the temperature (in K) obtained for the bin, and n is 
the corresponding total density (electron density plus proton den- 
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Figure 3. The Lx '■ Tx relation for the groups in our sample, using Chan- 
dra data. Temperatures and luminosities were extracted in regions that ex- 
tended out to 0.05i?5on- The solid points indicate the radio loud groups 
and the open points the radio quiet groups. We define radio loud groups as 
groups whose BGG has log(Li4on) > 21.5. In both panels, the solid line 
is the fit to the Lx '■ Tx relation for the entire sample (see Section l3~6t , the 
dotted line is the fit to the radio quiet sample and the dashed line is the fit to 
the radio loud sample. In the top panel, we fit gradients and intercepts for 
all three plotted relations, whereas in the bottom panel, we assume that the 
gradient of the relation for the entire sample holds for both the radio quiet 
and the radio loud sample, and fit only the intercepts. 



sity) for the bin. Here S is not a measure of the true entropy, but 
is rather the adiabatic constant - the ratio of the heat capacity for 
heating at constant pressure to the heat capacity for heating at con- 
stant volume (Cp/Cv), and is linked to the entropy of the gas (S r ) 
by 



S = 



h 2 



— exp (S'/cv - 5/3) , 



(4) 



27T (/it m H ) 

where h is the Hubble constant scaled to Ho = 
100 km s _1 Mpc -1 , fi is the reduced mass of hydrogen, 
tuh is the mass of a hydrogen atom, and c„ is the specific heat 
capac ity for expansion at constant volume (see iBalogh et al.l 
[19991 and references therein). Thus, S provides the most direct 
observational indicator of gas entropy. 

Pressure profiles P(r) are generated similarly using 

P(r) = [2.98 x 10 _10 Pa keV" 1 cm 3 ] n p (r)T(r), (5) 

In calculating the cooling time, T coo i, we note that as most of the 
groups have temperatures <2 keV , the thermal bremsstrahlung ap- 
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Figure 4. X-ray temperature profiles for the groups in the GEMS sample 
which have Chandra data. Groups containing a BGG with logZ/1400 > 
21.5 are plotted with closed points, and groups with BGGs with a radio 
power less then this are plotted with open points. 



Table 3. Results of power law fits to our scaled temperature, density and en- 
tropy profiles. The normalisations represent the logarithmic values for 
temperature, density and entropy of a 1 keV group at K500J assuming 
that the power law remains unbroken to this distance. 



Quantity 


Index 


Normalisation 


Temperature (keV) 


0.26±0.05 


0.36±0.05 


Density (cm -3 ) 


-1.23±0.06 


-4.1 ±0.1 


Entropy (keV cm 2 ) 


l.O8±0.05 


3.05±0.10 



NOTES: The sub-samples were fitted in logarithmic space using an or- 
thogonal regression method. Errors are 1 it for two interesting parame- 
ters. 



proximation is not valid. Instead, an estimate of the cooling time 
can be made by dividing the total internal energy of the gas in a 
given radial bin 



U(r) = ~(n e 



■n p )VT(r), 



(6) 



by the bolometric luminosity, L(, ;, obtained from XSPEC when fit- 
ting the data. 

However, since the groups are all of different masses, a mean- 
ingful comparison of the derived quantities cannot be made with- 
out scaling to take the mass di fferences into account . Here, we use 
modified self-similar scaling JVoit & Ponmanll2003l) such that ra- 
dial distances are scaled by R500 (the radius at which the density of 
the system is 500 times the critical density of the universe), temper- 
atures are scaled by T ave , the avera ge temperature of the syste m (as 
obtained from the ROSAT data of (Osmond & Ponmanll2004T) with 



1/2 



cool cores excluded from the spectral analysis), density by T a (, 

2/3 3/2 
entropy by T a „ e , pressure by T a „ e . We do not scale the cooling 

time profiles since in the regime where T <2 keV X-ray emission 

is dominated by line emission which is a function of temperature, 

which makes scaling by a single temperature unrealistic. The scaled 

profiles are shown in Figs l4l8l 

The scaled temperature, density and entropy profiles are then 

fitted with power laws using an orthog onal regression algorithm a s 

implemented in the SLOPES package dAkritas & Bers hadvlll99"r3) . 

The results of the fitting are shown in Table [3] 
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Figure 5. Gas density profiles for the same groups as in Fig[4] The data are 
scaled horizontally as in Fig|4]and vertically by T 1 ' 2 . 
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Figure 8. Gas cooling time profiles for the same groups as in Fig [4] The 
data are scaled horizontally as in Fig|4]and vertically unsealed. 
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Figure 6. Gas entropy profiles for the same groups as in Fig|4] The data are 
scaled horizontally as in Fig|4]and vertically by T 2 / 3 . 
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Figure 7. Gas pressure profiles for the same groups as in Fig [4] The data 
are scaled horizontally as in Fig|4]and vertically by T 3 / 2 . 



3 RESULTS 

3.1 General observations 

Examining the scaled profiles, it appears t hat groups follow mod - 
ified self-similar scaling, as discussed by [vbit & Pon manl d2003h . 
Looking closely at the entropy profiles, it is apparent that there is 
no isentropic core, even down to small fractions of the virial ra- 
dius; the entropy profiles follow predictions from models of cool- 
ing of gas without any significant heating to raise the central en- 
tropy. Thi s is in contrast to work probing hotter galaxy clusters, 
e.g. that of lVoit & Donahue] d2005t) who find isentropic cores within 
10 kpc(~ O.OI-R500) in 2.2-5.5 keV clusters. The profiles pre- 
sented in this paper probe the hot gas in groups to similar fractions 
of -R500, but no evidence of isentropic cores is seen. Sanderson et al 
(in prep) also find a similar result to that presented here for galaxy 
clusters observed with Chandra . 



3.2 Radio loud versus radio quiet sources 

The observed profiles are of a mixture of groups with and without 
current radio activity in the BGG. If radio outbursts are injecting 
energy into the IGM, then we might expect these outbursts to have 
an effect on the gas properties of the IGM. In order to investigate 
this, we divide our sample into two, based on the 1 .4 GHz radio lu- 
minosity (log [L1400 / (W Hz" 1 )] ) of the central AGN ( given in 
Tabled- We split our sample at three different radio luminosities, 
log (I/moo ) =21, 21.5, and 22, in order to establish if there are any 
differences as radio power increases. The scaled density, tempera- 
ture, and entropy profiles for each sample are then fitted with power 
laws as described in Section l2~3l 

Fitting the subsamples, we find that whilst there are some 
differences between the gradients of the scaled temperature and 
entropy profiles, the differences are slight, approximately l-2r. 
Whilst this could indicate a possible difference between the tem- 
perature and entropy profiles of radio loud and radio quiet groups, 
we cannot rule out that it is due to statistical scatter. To test this 
further, the profiles were co-added and overplotted to see if the av- 
eraged profiles showed any significant differences. The averaged 
profiles are shown in Fig[9]for the L1400 =21.5 cut, and the results 
of the fits for all three cuts are shown in Table [4] It can be seen that 
our data show that there is approximately a 2a difference between 
the temperature gradients of the radio loud sample compared to the 
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Table 4. Fits to our sample split according to the radio power of the BGG. 



Sub-sample Temperature Density Entropy 

gradient intercept gradient intercept gradient intercept 



^1400 


< 21 


0.22±0.07 


0.3±0.1 


-1.2±0.1 


-4.1 ±0.1 


l.ldbO.l 


3.0±0.2 


-£-1400 


> 21 


0.28±0.03 


0.38±0.05 


-1.23±0.07 


-4.1 ±0.1 


1.09±0.07 


3.1±0.1 


-£-1400 


< 21.5 


0.24±0.04 


0.38±0.06 


-1.16±0.07 


-3.9±0.1 


1.00±0.08 


2.9±0.2 


-£-1400 


> 21.5 


0.30±0.01 


0.42±0.02 


-1.08±0.09 


-3.7±0.2 


1.04±0.05 


3.02±0.09 


-£-1400 


< 22 


0.23±0.04 


0.33±0.07 


-1.17±0.07 


-3.9±0.1 


1.01±0.08 


2.9±0.1 


-£-1400 


> 22 


0.34±0.04 


0.45±0.06 


-1.3±0.1 


-4.2±0.2 


1.21±0.09 


3.2±0.1 



NOTES: The sub-samples were fitted in logarithmic space using an orthogonal regression method. Errors are 1<t for two interesting parameters. 



radio quiet sample (a gradient of 0.31±0.01 for the co-added radio 
loud sample compared to 0.24±0.04 for the co-added radio quiet 
sample). The apparent jump in temperature seen at O.OO4.R500 
in the radio quiet sample does not seem to be a feature associ- 
ated with all radio-quiet groups, but rather primarily the result 
of one anomalous system, namely NGC 3607. No significant dif- 
ference in gradients is seen between the co-added scaled radio loud 
and radio quiet entropy profiles (gradients of 1.04±0.05 for the ra- 
dio loud sample compared to 1.00±0.08 for the radio quiet sam- 
ple). Further, the radio quiet and radio loud density profiles show 
little difference in gradient ( — 1.16 ± 0.07 for the radio quiet com- 
pared to —1.09 ± 0.09 for the radio loud samples). 



3.3 Black hole mass and heating 

As only a small difference is seen between radio quiet and radio 
loud systems, it could be that the effect of an AGN is cumulative, 
building up over repeated duty cycles. If this is so, then comparing 
profiles of groups with more massive central black holes with those 
groups with smaller black holes could provide some insight, since 
a larger black hole should theoretically have had more accretion 
cycles than a smaller one, and thus, more cycles of radio activity 
that should have injected more energy into the group. 

The relationship between black hole mass and galactic 
bulge mass has been investigated by many authors (for example 
M arcom & Huntl2003l) . and whilst such a relationship does exist, it 
contains significant amounts of scatter, thus we must exercise cau- 
tion in interpreting correlations that make use of it. 

We obtain K-band luminosities (Mk) for the BGGs in our 
sample from 2MASS (see also Table[2}, and using Mk as a proxy 
for black hole mass, we divide our sample into two based on the 
Mk of the BGG, and create subsamples for Mk < —25 and 
Mk > —25. We co-add the scaled temperature, density, and en- 
tropy profiles for the subsamples, and fit them as in Section \32\ 
The results of fitting the scaled profiles of the two subsamples are 
shown in Table[5] and the co-added profiles are presented in FigllOl 

From the plots it can be seen that there is a difference in the 
scaled temperature profiles in that groups which host larger BGGs 
(hence those BGGs which should have more massive black holes 
and have undergone more cycles of accretion) have a lower scaled 
temperature within 0.047?soo than groups with smaller BGGs. This 
could arise from both subsamples having the same central temper- 
ature but with the Lk bright subsample having an higher over- 
all group temperature. Fitting the binned profiles, we find that the 
temperature gradient for the sample containing groups with larger 
BGGs is steeper than the sample of those groups with smaller 



BGGS (0.28±0.03 compared with 0.22±0.02). From Table|5]it can 
be seem that there is a small (approximately 1. 5a) difference be- 
tween the entropy profiles of the two subsamples; the fits to the 
binned entropy profiles show approximately a la difference be- 
tween the slopes of the entropy profiles (1.04±0.07 for the sample 
withMx < -25 and 0.96±0.08 for the sample with M K > -25). 
This could be a result of a simple correlation between the X-ray 
scaled gradients of the hot gas and the BGG mass, since larger 
BGGs would have steeper underlying potential wells which would 
lead to steeper temperature gradients if the X-ray gas is a tracer of 
the dark matter. Further, we cannot rule out that any correlation be- 
tween our scaled temperature, density and entropy profiles may be 
a result of the size of the host galaxy rather than black hole mass. 



3.4 Using an alternative scaling 

In Sections [3.2l and l33l the gas profiles have been scaled by an av- 
erage temperature obtained from the ROSAT data. As we expect the 
Chandra temperature profiles to agree with the ROSAT profiles at 
large distances from the centre of the group, scaling the Chandra 
profiles using the ROSAT temperature almost guarantees conver- 
gence of scaled profiles outside the cores of groups. An indepen- 
dent temperature scaling would be desirable, to allow an absolute 
normalisation of profiles and avoid introducing an artificial conver- 
gence. To this end, we investigate the use of the velocity dispersion 
of the galaxy group, ac as a proxy of the virial temperature Ty. 
Since TV oc ffg, scaling the entropy profiles by er^/ 3 should be 

2/3 

equivalent to scaling by T v . 

We can also improve on the use of Mk as a proxy for the 
mass of the central black hole. The Mk -black hole mass (Mbh) 
relation shows substantial scatter, whilst the a : Mbh rela- 
tion, where a is the velo city dispersion of the BGG, is tighter 
jFerrarese & Merrittl20 00 and Gebhar dTet alj2000h . However, re- 
cent work (Lauer et al. 2006) suggests that for BGGs, the a : 
Mbh relation may not be as well constrained compared to the 
My : Mbh relation, where My is the V-band magnitude of 
the BGG. This is however a much debated issue, and it has 
been suggested that the discrepancy between Mbh as predicted 
by the two different methods a rises du e to My being overesti- 
mated in BGGs (Batcheldor et al. 2006). They argue that in fact 
when black hole masses obtained using a a : Mbh relation are 
compared to those using a Mk ■ a relation, i.e. a n ear-infra- 
red r ather than blue magnitude relation, (e.g. [Marconi & Hunt! 
2003), there is no discrepancy, and that the discrepancy arises 
due to the unusual colour profiles of BGGs compared to E and 
SO galaxies. With this in mind, and with no apparent resolution, 
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Table 5. Fits to our sample split according to the mass of the central black hole. 



Sub-sample 


Temperature 


Density 


Entropy 




gradient intercept 


gradient intercept 


gradient intercept 


M K < -25 


0.28±0.03 0.43±0.07 


-1.30±0.07 -4.1 ±0.1 


1.04±0.07 3.1 ±0.1 


M K > -25 


0.22±0.02 0.36±0.06 


-1.30±0.08 -4.1 ±0.1 


0.96±0.08 3.1=1=0.1 



NOTES: The sub-samples were fitted in logarithmic space using an orthogonal regression method. Errors are 1<t for two interesting parameters. 



we opt to use the tighter a : Mbh relation, and scale the entropy 
profiles by Rsoo radially, and by a^ 3 vertically, and split the sam- 
ple into two in two different ways; in Figll lh. as in Section [3~2l at 
log (Li4oo) = 21.5 and in Figll lb at a = 270 km s _1 . Values for 
a are collated from Hyperleda_j, and shown in Table [2] The sub- 
samples for both the log (L1400) and a splits are fitted as before. 
We find that the entropy profiles in the log (L1400) split show no 
significant differences between the profiles, whilst the a split shows 
that systems with higher values of a have steeper entropy profiles 
than those with lower values (gradients of 1.22±0.0T compared to 
0.78±0.05), suggesting that the effect seen in Section [3~3l is real. 

3.5 Effects of AGN activity on individual systems 

As has been discussed in Sections [3.2l and [331 when comparing the 
sample as a whole, there appears to be some marginal difference 
between the gas properties of the groups with radio loud BGGs 
and those with radio quiet BGGs. The fact that the difference is 
marginal could indicate that a real effect is being diluted by aver- 
aging together groups with different radio properties. To test this 
we correlate the logarithmic gradients of scaled temperature, en- 
tropy and density profiles of individual groups with log(Li4oo) 
of the BGG The Pearson R correlation coefficients are shown in 
Table [6] and plots of the various quantities that have been corre- 
lated are shown in Fig.[T2] We find significant correlations between 
L1400 and both the temperature and entropy gradients. 

Since our sample contains three systems with extensive jet- 
fed radio lobes, it is interesting to examine whether these systems 
might be entirely responsible for the rather weak correlations we 
see. These three systems (NGC 383, NGC 741, NGC 4261) have 
jets which transport the AGN energy out to large radii where the 
gas density may be too low to allow the energy to be radiated in 
less than a Hubble time. This could lead to steepening of tempera- 
ture and entropy profiles in these systems. We find that if the three 
groups with large radio galaxies are excluded from the sample the 
correlations are no longer significant. Hence our data do not re- 
quire any effect on the gas from AGN beyond that arising in the 
three large jet systems. 

We further correlate L1400 of the BGG, with the global X-ray 
temperature, Tx, as measured by ROSAT , and with Mk, the K- 
band magnitude of the BGG (Table [33] and Fig.|13t. We also find 
that L1400 correlates with Tx- This correlation suggests that either 
the radio source does heat the IGM, or alternatively, that a hotter 
group, with a strong temperature and entropy gradient is more con- 
ducive to triggering an AGN outflow of the scale of the largest ra- 
dio galaxies, than a cooler group. There appears to be no significant 
correlation between L1400 and Mk or Mk and Tx, indicating that 
the effect is not due to larger galaxies hosting larger radio outbursts 
and being found in hotter groups. This suggests that the difference 

1 http://leda.univ-lyonl.fr 



Table 6. Correlation coefficients for logarithmic temperature, entropy and 
density gradients and log (I/i40o)- The final column indicates whether the 
correlation is significant at the 95% confidence level, i.e. can the null hy- 
pothesis (that there is no correlation) be rejected with only a 5% level of 
uncertainty. 



Quantity 




R 


Significant? 


dlogT/dlof 


;R 


0.50 


Y 1 


d log S/d log 


;R 


0.56 


Y 1 


d log n/d log 


;R 


-0.23 


N 



NOTES: x The correlation disappears when the sources with powerful 
radio jets are removed from the sample (see text). 

Table 7. Correlation coefficients for radio luminosity, X-ray temperature, 
and Mg of the systems in our sample. Correlations marked with 1 are those 
for which a correlation significant at 95% was found. 

Quantity L1400 T x M K 

L1400 0.51 1 -0.14 

T x - - 0.46 



in temperature profiles found in Section [3~3l mav arise from steeper 
potential wells in the systems with larger BGGs. 

3.6 The L x : T x relation 

ICroston et al] (2005) find that groups which host radio galax- 
ies appear to fall below t he Lx ■ Tx relation for groups (c.f. 
lOsmond & Ponmar] |2004;)), which they interpret as evidence for 
AGN heating of galaxy groups. We use the temperatures and lumi- 
nosities extracted in Section |2.3.1| to produce an Lx ■ Tx relation 
for the inner regions of groups as sampled by the Chandra data; the 
corresponding Lx '■ Tx relation is shown in Fig. [3^. 

From Fig. [3] there does appear to be some difference in the 
distribution of radio loud and radio quiet points around the mean 
Lx '■ Tx relation, and to test this, we measure the perpendicu- 
lar distance for each group from the mean Lx ■ Tx fit, and as 
above, compare the distribution of distances for the radio loud and 
radio quiet systems. A K-S test indicates that there is a 20% chance 
that the observed scatter around the mean relation would arise by 
chance if the radio loud and radio quiet sub-samples came from 
the same population (D=0.46, P=0.20), implying that the effect is 
marginally significant. However, we are dealing with only 19 sys- 
tems, and so both effects need to be explored with a larger sample. 

If the difference in distributions around the mean relation is 
real, then fixing the gradient of the relation to that obtained from the 
fit to the entire sample allows an estimate of the temperature offset 
between the radio loud sample and the radio quiet sample to be de- 
termined. We find that the temperature offset AT = 0.2 ± 0. 1 keV. 
If this temperature offset is entirely due to the impact of AGN heat- 
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Figure 9. Co-added, self-similarly scaled temperature, density and entropy 
profiles (from top to bottom), of the sample split according to the radio 
power of the BGG. The closed points indicate groups with BGGs such that 
£l400 >21.5, and the open points comprise the rest of the sample. There 
is very little difference between the density and entropy profiles of groups 
containing radio loud BGGs and those containing radio quiet BGGs, but 
that there is a significant (2a) difference between the gradients of the tem- 
perature profiles. 



ing within 0.05i?soo, then an estimate of the energy required to 
produce this offset can be obtained. 

The energy input AE required to produce a temperature 
change of AT is given by: 



AE = NkAT 



(7) 



Figure 10. Co-added, scaled temperature, density and entropy profiles 
(from top to bottom) of the sample split into two subsamples, based on Afjf . 
The solid points represent the subsample for which Mk <-25, whereas the 
open points represent the subsample where Mk >-25. There appears to be 
a difference between the temperature profiles of the two subsamples. 



where TV is the total number of particles within 0.05-R 5 oo and k 
is the Boltzmann constant. Assuming that the average radio quiet 
profile is roughly equivalent to an unheated group, we use the den- 
sity profile shown in Fig|9p, and a AT = 0.2 keV to estimate AE 
for a 1 keV system. We find that AE ~ 2.6 x 10 57 erg. Assuming 
an AGN lifetime of 10 7 yrs, this gives an energy injection rate of 
~ 8 x 10 42 erg s~ x . We comment below on the plausibility of this 
injection rate. 

We further compare the marginal effect found in the Chandra 
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Figure 11. Entropy profiles for the sample split and binned according to 
log (L1400) (top) and the galaxy velocity dispersion, a, (bottom). In the 
top plot, the open points represent those groups in which log (L1400) < 
21.5 and the black points are those where the BGG radio power is 
greater than this limit. In the bottom plot, the open points represent 
those systems where a < 270 km s -1 , and the black points are the 
systems where a > 270 kms -1 . The profiles are scaled radially by 
-R500 and vertically by cfq 3 . We see very little difference between the en- 
tropy profiles of radio loud versus radio quiet groups, whilst groups whose 
BGGs have large black holes have steeper entropy profiles than those groups 
whose BGGs have smaller black holes. 



data with the result of ICroston et al. (2005), to ascertain whether 
any AGN heating is taking place within the core, compared to out- 
side the core. To do this, we examine the mean residuals (<tT) of 
the radio loud and radio quiet subsamples co mpared to the fit to th e 
radio quiet sample. For the ROSAT sample o flCroston et alj J2005h . 
we use the best fitting least squares fit to their c2 radio quiet sample 
in order to have a reasonable comparison of radio power with our 
log (Li4oo ) = 21.5 cut. We find that in the ROSAT data, the mean 
residual for the radio quiet group is <jTi?Q=0.06±0. 1 keV, com- 
pared to the value for the radio loud sample, gTrl=1-6±Q.5 keV. 
The corresponding quantities (ET) for the Chandra data are 
T,T RQ =0.07±0. 1 keV and TT RL =0.2±0. 1 keV. The mean residual 
is much larger for the radio loud sample in the ROSAT data than in 
the Chandra data, implying that more significant AGN heating may 
be taking place outside the region being probed here (~ O.I-Rsqo). 
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Figure 12. The logarithmic density, temperature and entropy gradients plot- 
ted against Li4oo- There exists some correlation between the radio power 
and the entropy and temperature gradients, but not with the density gradient. 



4 DISCUSSION 

To summarise the key results above, we find modest differences 
between the gas properties of the cores of radio loud groups com- 
pared to radio quiet groups; radio loud groups appear to have 
slightly steeper temperature profiles, and there are marginally sig- 
nificant correlations between the temperature and entropy gradient 
of groups and the radio power of their BGGs, though this may arise 
from the three strong jet sources in our sample. We also find that 
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Figure 13. Plots of Tx against L1400 (top), Mk 
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related. 



when the sample is split according by a proxy of the black hole 
mass (either Mk or a), steeper temperature and entropy gradients 
are found for those groups whose central galaxies have large values 
of Mk or a, compared to those groups where Mk or a are lower. 
Considering the effect of active radio galaxies on the scaling rela- 
tions, we find evidence for a small offset in the Lx '■ Tx relation 
between radio loud and radio quiet groups, b ut this is much smalle r 
than the offset seen on larger spatial scales bv lCroston et al.lf2005t) . 



In light of these results, we now discuss the three possible 
explanations for any observed correlations that were highlighted in 
SectionQ] 

• The radio source modifies the X-ray properties through 
heating or displacement of the gas. From the images presented 
in Section [2~2l and the work of other authors, it is apparent that 
the radio source can displace gas, creating cavities in the X-ray 
emission. The entropy profiles when split on radio power imply 
that no extra energy is being injected into the core by current ra- 
dio source heating, although systems which exhibit powerful jets 
may be transporting extra energy to regions outside the core, and 
causing a steepening of the entropy gradient outside the core. How- 
ever, splitting the sample using proxies of black hole mass suggests 
that there may be a cumulative effect of repeated AGN outbursts. 
It is not clear whether the cumulative effect is from repeated AGN 
outbursts, or due to energy injection from the stellar content of the 
galaxy. 

• The observed X-ray properties provide suitable conditions 
for radio source triggering. If this were the case, then it would 
be expected that the immediate gaseous environments of currently 
radio active sources would be different from the environments of 
currently radio quiet sources. It can be seen from Fig [9] that the 
profiles of radio loud and radio quiet groups are rather similar, and 
from Figs 141 8 1 the scatter between individual profiles is larger than 
any systematic differences seen in the averaged profiles. As any 
structural differences conducive to the production of a radio source 
would have to be most significant in the core of the group, since 
that is where the black hole resides, the similar profiles for the radio 
loud and radio quiet groups argue against this hypothesis. 

• There exists some third parameter that correlates with 
both radio source activity and position on the Lx ■ Tx rela- 
tion. From Fig 1 13b. it could be the case that Mk may play some 
role, as the mean magnitude for the radio loud sample does seem to 
be larger than the mean magnitude for the radio quiet sample. This 
could be a result of a cuspier potential caused by a larger central 
galaxy. 

The evidence suggests that our first and third hypotheses 
could, together, help explain the correlations found here, whilst 
the second hypothesis seems unlikely, since any systematic differ- 
ences between radio loud and radio quiet groups found are small, 
and can be attributed to either statistical scatter, or to a small num- 
ber of sources with physically large sizes that transport energy out 
to larger radii. The fact that there are no strong differences be- 
tween radio-loud and radio-quiet groups implies that in general ra- 
dio sources do not inject enough energy into the central regions of 
their host groups to cause observable heating effects. If there is any 
heating going on on these scales, it must just be enough to balance 
cooling. 

To investigate whether AGN output is indeed sufficient to 
balance cooling, we consider the case of NGC 383 (3C31). This 
source has been well studied and is one of the few in our sam- 
ple that has an estimated total jet power. We define the cooling ra- 
dius as the radius at which the cooling time is equal to the Hub- 
ble time, which for NGC 383 is ~ 25 kpc. Within this radius, 
Lboi ~ 5 x 10 41 e rg s -1 . The total jet e nergy flux at 12 kpc is 
- 9 x 10 43 erg s" 1 jLaing & Bridiel2 002). Further, in Section lTol 
we calculated that the average temperature offset between the ra- 
dio loud sample and the radio quiet sample is ~0.2 keV. Over 10 7 
years, an energy injection rate of ~ 8 x 10 42 erg s _1 is required to 
produce the offset in the NGC 383 group. This is about an order of 
magnitude less than the maximum AGN energy input rate, so that 
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even a less luminous AGN, where the AGN energy flux would not 
be as high, could provide sufficient energy to both counteract ra- 
diative cooling, and raise th e temperature of the gr oup. This agrees 
well with the conclusions of Ha rdcastle et al.l J2002I) . 

Further, the large scale radio morphology of NGC 383 and 
the fact that the bolometric luminosity of the cool core is much 
less than the radio source power implies that an outburst such as 
that in NGC 383 could transport much of the energy contained in 
the radio source out to larger radii, where i t could also affect the 
hot gas, explaining the lCroston et alj J2005!) results. However, the 
radio morphology of NGC 383 is very different to most of the other 
radio outbursts in this sample, and its long jets and plumes may be 
transporting energy to regions outside the core of the group without 
affecting the innermost regions. Radio outbursts on smaller scales 
may affect the IGM closer to the BGG. 



5 CONCLUSIONS 

We have presented azimuthally averaged profiles of the properties 
of the central hot gas, and an Lx ■ Tx relation within 0.05i?soo 
for 15 galaxy groups with varying degrees of AGN activity. 

The slightly steeper temperature profiles seen in radio loud 
groups could be related to the size of the BGG - a larger BGG 
would have a steeper potential well associated with it, and hence a 
steeper gradient - or it could be that systems with radio-loud AGN 
have cooled more than those without radio-loud AGN, and it is this 
cooling gas that fuels the radio source. However, our sample size 
is small, and the properties of the gas and how they correlate with 
the AGN activity need to be investigated on a wider scale using 
a larger sample which probes a wider range of temperatures and 
scales of AGN activity. Alternatively, the effect could arise from 
the contribution of a few radio sources with large physical sizes. 
These sources may have little effect on the gas close to the core of 
the group, depositing their energy at larger radii, at the outside edge 
of the regions probed here, causing the gas profiles of their host 
groups to appear to be marginally steeper over our radial range. 

The steepening of the entropy gradients in systems with larger 
values of Mk and a (indicating larger black holes) compared to 
those with smaller values of Mk and a could be evidence for re- 
peated cycles of AGN activity having an impact on the hot gas in 
groups. However, it is unclear whether the effect has to have oc- 
curred via AGN heating, rather than via other forms of energy injec- 
tion from the galaxy. As both AGN and galaxy growth are linked, 
disentangling one effect from the other may be difficult, particu- 
larly as AGN life-cycles are also poorly understood. Work needs 
to be done here in both comparing the gas properties of groups at 
different stages of their evolution, and on the life-cycles of AGN 
hosted by BGGs. 

The Lx '■ Tx relations show that any effect that AGN are hav- 
ing on the gas primarily occurs at larger radii; it could be that within 
the region probed by Chandra , AGN only act to counteract radia- 
tive cooling. The difference between the radio loud a nd radio quiet 
subsam ples is smaller than the difference found by ICroston et al.l 
j2005h . suggesting that ongoing large scale outbursts such as those 
in NGC 383 and NGC 4261 could be transferring energy out to 
large radii, and may have a stronger effect outside the core, whilst 
the effects of repeated outbursts are more likely to show up in the 
core. However, more X-ray work is needed here to further probe 
and constrain the gas properties at larger radii. 

We conclude that it appears that radio loud AGN do not ir- 
reversibly raise the entropy in the core of their host galaxy groups. 



Rather, it seems that the feedback mechanisms at work result in typ- 
ical central AGN generating just sufficient energy to balance cool- 
ing in the core of the group. Repeated outbursts may have some 
longer lasting effect on the gas, which may contribute to similar- 
ity breaking, although detailed studies extending to larger radii are 
required to investigate this. This cumulative effect may take place 
primarily via rather rare large outbursts, which dump energy out- 
side the cooling radius via radio jets. 
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